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Abstract 

We present robust algorithms for set operations and Euclidean trans- 
formations of curved shapes in the plane using approximate geometric 
primitives. We use a refinement algorithm to ensure consistency. Its 
computational complexity is 0(nlog?i + k) for an input of size n with 
k — 0(n 2 ) consistency violations. The output is as accurate as the geo- 
metric primitives. We validate our algorithms in floating point using se- 
quences of six set operations and Euclidean transforms on shapes bounded 
by curves of algebraic degree 1 to 6. We test generic and degenerate in- 
puts. 

1 Introduction 

Set operations and Euclidean transformations are important computational ge- 
ometry tasks with many applications. There are efficient algorithms for planar 
shapes, which are more common and simpler than 3D shapes. A shape is mod- 
eled as a subdivision: a partition of the plane into faces, curves, and points. 
Fig. QJi shows a subdivision with faces /o and /i, curves ei, . . . , e^, and points 
indicated with dots. A set operation is performed by constructing the mutual 
refinement of the input subdivisions, called the overlay, and returning the faces 
that satisfy the set operator (Fig.QJ). A Euclidean transformation is performed 
by transforming the points, curves, and faces of the input subdivision. 

Shape manipulation algorithms use geometric primitives that are formulated 
in the real- RAM model where real arithmetic is exact and has unit cost. The ro- 
bustness problem is how to implement these primitives in computer arithmetic. 
The mainstream strategy is to implement primitives exactly using integer arith- 
metic and algebraic computation. We prefer approximate primitives in floating 
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(a) (b) (c) 

Figure 1: Subdivisions (a, b) and their overlay (c) with fi n /2 shaded. 

point because they are much faster and have constant bit complexity. Although 
the approximation error is negligible, it can cause prior shape manipulation 
algorithms to generate inconsistent or highly inaccurate outputs. We present 
shape manipulation algorithms that use approximate primitives yet generate 
consistent, accurate outputs. 

1.1 Prior work 

Yap[JJ describes exact computational geometry. Exact computation increases 
bit complexity, hence running time. The CGAL library[5] provides floating 
point filtering techniques that somewhat reduce this cost. The increase in bit 
complexity from input to output leads to unbounded complexity in sequences 
of computations. The only current solution is to simplify the output. Although 
simplification algorithms have been developed for a few domains, [U EJ 16] not 
including implicit algebraic curves, a general strategy is unknown. 

The alternative to exact primitives is approximate primitives. The robust- 
ness problem is to ensure consistency despite approximation error. 

The controlled perturbation robustness strategy is to verify that the pred- 
icates in an algorithm are correct despite approximation error. Verification 
fails when the predicate value is too close to zero. Failure is prevented prob- 
abilistically by randomly perturbing the input parameters. If any predicate 
still cannot be verified, the input is perturbed again. Controlled perturbation 
has been applied to arrangements, [9[ QUI [11] convex hulls, [12] and Dclaunay 
triangulation. [12] The perturbation size is exponential in the algebraic degree 
of the primitives. The degree-four perturbation size exceeds the error bounds of 
engineering applications. We developed versions of controlled perturbation that 
solve this problem and used them to compute Minkowski sums of polvhedra|13| 
and free spaces of planar robots. [14] Controlled perturbation does not work for 
shape manipulation because it does not handle algebraic curves and does not 
address output simplification in sequences of computations. 

Our robustness strategy is to eliminate inconsistency due to approximate ge- 
ometric primitives by enforcing consistency constraints on the data structures 
that represent geometric objects. We analyze the running time and the accuracy 
in terms of the approximation error, S, and the number of constraint violations, 
k, but neither parameter is provided to the algorithms. An algorithm is incon- 
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sistency sensitive when the extra running time for inconsistency elimination is 
polynomial in k and n, the extra error is polynomial in k and n and is linear 
in 6, and k is polynomial in n. Inconsistency sensitivity captures the concept 
of an efficient approximate algorithm: fast and accurate when there are few 
inconsistencies, and degrading gracefully 

We developed an inconsistency sensitive arrangement algorithm for algebraic- 
curve segments. [7] The x-coordinates where segments cross and the segment y- 
order between crossings are computed with (1 + kn)e accuracy for n segments, 
where e denotes the error in the approximate geometric primitives of that algo- 
rithm. A consistency constraint violation is a cyclic y-order, due to an incorrect 
ordering of segment crossing coordinates. The algorithm computes a consistent 
y-order of size V = 2n + N + mm(3kn, n 2 /2) for TV crossings. The running time 
is 0(V\ogn) and the output is correct for a set of segments within (1 + kn)e of 
the output segments. We also developed two inconsistency sensitive Minkowski 
sum algorithms. [5] 

The arrangement algorithm can compute a curve y-order that is inconsistent 
with the curve endpoints, as illustrated below. In prior work, we developed a 
simple solution for the special case where every curve endpoint has a distinct 
x-coordinate. We handled the general case heuristically by sweeping along a 
random axis. This strategy proved inadequate even for generic inputs. In this 
paper, we present a general, inconsistency sensitive solution. 

1.2 Inconsistency in the Overlay Algorithm 

Our overlay algorithm uses three geometric primitives: computing intersection 
points of monotone curves, computing turning points of non-monotone curves, 
and comparing coordinates of points. The first primitive is required by any over- 
lay algorithm, while the other two are dictated by our use of a sweep algorithm. 
The high level algorithm is as follows. Split the curves at their turning points. 
Split the monotone curves at their intersection points to obtain sub-curves. De- 
rive the y-order of each pair of sub-curves that overlap in x by comparing the 
y coordinates of their intersection points with a vertical line at the middle of 
their interval of x-ovcrlap. The details of the sweep and how it avoids creating 
cycles appears in previous work[5] and in Sec. [5] 

The algorithm output, the x-order of the curve endpoints and the partial y- 
order of the sub-curves, allows us to compute the subdivision structure without 
additional approximate primitives. In particular, we can use a standard sweep 
algorithm to generate a vertical cell decomposition. The curve endpoints are the 
sweep events in order of increasing x. The sweep list contains curves in y-order 
whose left endpoint event has been processed and whose right endpoint has not. 
If a subdivision construction algorithm, e.g. randomized incremental vertical 
cell decomposition, requires endpoint/curve order, the y-order of endpoint a of 
curve ab with respect to curve cd, with c x < a x < d x , is deduced from the 
calculated y-order of ab and cd. 

We illustrate that an inconsistent y-order can arise in computing the overlay 
of a monotone curve, e2, with two monotone curves, e\ and e3, that meet at a 
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Figure 2: Inconsistent y-order (a), e2 deformation (b), and consistent refinement 
(c). 

point, b (Fig. Due to numerical error, the intersection primitive incorrectly 
reports that e-i and do not intersect. Intersections with vertical lines (arrows) 
order ej above e-i and above e%. The incorrect y-order of e-i and is correct 
for a small deformation of the input (Fig. ^jp). But the algorithm output is 
inconsistent because no deformation can make e-i be below e\, above e^, and 
disjoint from b. 

1.3 Contribution and Organization 

Our main result is an inconsistency sensitive refinement algorithm that assigns 
a consistent y-order to a set of curves by splitting some curves at endpoints 
of other curves (Fig. [2]:). The output y-order is the induced refinement of the 
input y-order. In our example, ei is split into e2i and ei2 at b, e2i is below ei, 
e22 is above and the other pairs are unordered because their domains are 
disjoint. The precondition is that each pair of input curves can be placed in 
the input y-order using a deformation of at most 8. The postcondition is that 
a single (5-deformation can place all the output curves in the output y-order. 
The running time is 0(nlogn + k) for n curves with k = 0{n 2 ) consistency 
violations. We define concepts in Sec. [2j present the refinement algorithm in 
Sec. [3J and analyze it in Sees. HHSj 

As indicated above, overlay can introduce inconsistencies. Surprisingly, even 
transformations can do so. To accomplish set operations, we start with a possi- 
bly inconsistent overlay, run the refinement algorithm, construct the faces, and 
extract the relevant faces (Sec. EJ. To transform a subdivision, we transform 
its curves, split rotated curves at new turning points, run refinement, and con- 
struct the faces (Sec. [7]). In this manner, the refinement algorithm enables us 
to accomplish consistent and accurate shape manipulation using approximate 
geometric primitives. 

We present an empirical validation that our algorithms arc fast and accurate 
on sequences of six set operations and Euclidean transformations on shapes 
bounded by curves of algebraic degree 1 to 6 (Sec. [5]). Moreover, k is zero for 
generic input and is small for degenerate input. We conclude with a discussion 
(Sec. Hi. 
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2 Concepts 



In this section, we define approximate subdivisions and four ways they can 
be inconsistent. Fig. [2^, is an example of an approximate subdivision and the 
depicted inconsistency is an instance of one of these ways. The refinement 
algorithm eliminates these inconsistencies. We prove that the output represents 
a consistent shape (Sec. 2]) and that its error is bounded by the input error 
(Sec. [SJ. Each group of related concepts is defined informally then formally. 

We work in the xy plane. A point, a, has Cartesian coordinates (a x ,a y ). 
Curves are open, bounded, and piecewise algebraic. They are assumed x- 
monotone and y-monotone unless stated otherwise. A curve, e, whose endpoints 
are a and b with a x < b x is denoted e[a, b]. It is vertical when a x = b Xl horizontal 
when a y ~ by, increasing when a y < by, and decreasing when a y > b y . 

We represent a subdivision as a set of disjoint curves (Dcf . [1} . The combi- 
natorial structure of the subdivision is derivable from the x-ordcr of the curve 
endpoints and the < y partial order on the curves (Dcf. [3]), as explained in the 
introduction. 

Definition 1 (subdivision). A subdivision is a set of curves in which every 
member is disjoint from the closure of every other member. 

Definition 2. Sets s and t overlap in x if there exist p € s and q S t with 
Px = q x - 

Definition 3 (< y ). For two point sets, s and t, that overlap in x, s < y t if 
Py < 1y f or eac/i p£s and q £ t with p x = q x . For a point, a, a < y s means 
{a} < y s. 

Wc model a shape with an approximate subdivision: a set of curves with an 
imposed acyclic partial y-order (Def. 2]) denoted -C y . An approximate subdi- 
vision is a subdivision when its curves are disjoint from the closures of other 
curves and the < y order matches the < y order. 

Definition 4 (approximate subdivision). An approximate subdivision is a set 
of curves with a binary relation, denoted -< y , whose transitive closure is irrefiex- 
ive, such that e -< y f or f -< y e iff e and f overlap in x. The notation e < y f 
means e < y f or e = f , e )~ y f means f -< y e, and e > y f means f -< y e. 

Fig. [3^, shows an example in which the -< y order for curves that overlap in 
x is index order, e.g. £3 ^ y e$ because 3 < 5. Although -< y agrees with < y 
here, they disagree elsewhere. For example, e^ intersects e^, while eg < y eio yet 
eio <y eg. 

An approximate subdivision (Fig. [3Ji) has a realization when an endpoint- 
preserving deformation (Fig. \Sjp) of the curves can make the actual < y order 
agree with the imposed < y order (Def. [5]). A realizable approximate subdi- 
vision represents an infinite number of shapes, but their curves all have the 
same endpoints and < y order, hence all these shapes have the same combinato- 
rial structure. An unrealizable approximate subdivision does not represent any 
shape. 
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(a) (b) 
Figure 3: Approximate subdivision (a) and realization (b). 



Definition 5 (realization). A realization of an approximate subdivision, S , is 
a function, r, from S onto a subdivision, T, such that for all e G S , e and r(e) 
have the same endpoints and for all e, / G S , e -< y f implies r(e) < y r(f). 

As indicated in the introduction, the vertical order of an endpoint, p, of a 
curve, /, with respect to a curve e, can be deduced from the order of e and /. 
In this manner, we define the lower and upper sets of endpoints, L(e) and U (e), 
with respect to e (Def.O, from the imposed vertical order -< y . We use these sets 
to define four consistency constraints (Def. [7]). Our central theorem is that an 
approximate subdivision has a realization iff it satisfies these four consistency 
constraints (Thm. [3]) . 

Definition 6 (Fie) and U{e)). For a curve, e[a, b] G S, an approximate sub- 
division, L(e) is the set of endpoints, p, of curves, f -< y e, with a x < p x < b x ; 
likewise, U(e) with f > y e. 

Definition 7 (consistent). An approximate subdivision, S, is consistent if each 
curve, e[a, b] £ S, satisfies four constraints. 

1. For each q G U(e), there is no p G L(e) with a x < p x = q x < b x and 
P V > q y . 

2. For each f[c,d] >- y e, a x = c x implies a y < c y and b x = d x implies 

by < dy. 

3. a) For each p G L(e), there are no c,d G U(e) with c x < p x < d x and 
c y , dy < p y . b) For each q G U(e), there are no c,d G He) with c x < q x < 
d x and c y ,d y > q y . 

4-. If e is vertical, there is no endpoint, p, of curve f , with p x — a x and 
p y G (a y ,b y ), the open interval between a y and b y . 

Fig. |4] illustrates the constraints. Part (a) violates constraint 1 if e -< v g, 
f -< y e, p x = q x , and p y > q y . (Fig. [2ji is a special case with e = en, / = e.3, 
g = ei and p = q = b.) Part (b) violates constraint 2 if e -< y /, a x = c x , and 
a y > Cy. Part (c) violates constraint 3a if f < y e, e < y g, e -< y h, c x < p x < d x , 
and p y > c y ,d y . Part (d) violates constraint 4 if and a y < p y < b y . 

All four of these violations occur in practice (Sec. [5]). 
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Figure 4: Inconsistent approximate subdivisions. 




Figure 5: Refinements of the approximate subdivisions in Fig. |4j 



Theorem 3. An approximate subdivision is realizable iff it is consistent (p. \17\ ). 

A refinement of an approximate subdivision is an approximate subdivision 
whose curves are sub-curves of the original curves with the induced -< y order 
(Def-EJ. 

Definition 8 (refinement of curve). A refinement of a curve, e[a, b], is a set of 

curves, C(e), that form an x-monotone chain from a to b. 

Definition 9 (refinement of approximate subdivision). A refinement of an ap- 
proximate subdivision, S, is an approximate subdivision, S' , with the same set 
of endpoints in which each curve, e G S, maps to a refinement, C(e) C S' , and 
e -< y f in S implies e 1 < y f in S' for every e' G C(e) and f G C(/) that 
overlap in x. 

Fig. [5] shows consistent refinements of the inconsistent approximate subdivi- 
sions from Fig.@]witli C(e) = {e\, e2, e^} in parts (a) and (c) and C(e) = {ei, e-{\ 
in parts (b) and (d). The refinements are realizations, except in part (c) where 
ei < y f but / -< y e, and h < y e% but e -< y h. However, a realization exists by 
Thm.[l 

Theorem 1. The refinement algorithm computes a consistent refinement fv. \l!i!\) . 

The remaining concepts relate to error analysis. We define an error metric 
on approximate subdivisions that models errors in the pairwise -< y orders due 
to errors in the geometric primitives used to compute them, especially curve 
intersection point computation. An approximate subdivision is 8-accurate when 
the -< y order of each pair of curves equals the < y order of a pair of <5-close curves 
(Def. HT|). Each pair of curves can be placed in -< y order by a small deformation. 

Definition 10. A ^-deformation of a curve, e, is an x-monotone curve with 
the same endpoints as e whose Hausdorff distance from e is less than 5. 
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Figure 6: Illustrations of 5-accuracy. 



Definition 11 (8- accurate). An approximate subdivision is 5-accurate if e -< v f 
implies that e and f have S- deformations, e and f , with e < y f . 

Fig. [B^-b show that Fig. 0k is 5-accurate for a small 6, which matches our 
intuition that the curves barely intersect. The pairs e, f and e,g use different 
e curves. Fig. \§jp shows that Fig. HJd is <5-accurate. Even though c <G e, e < y f 
because curves are open. Fig. [6jl shows that Fig. is (^-accurate. The curves 
are their own 5-perturbations, and e is x-monotone but not monotone. Fig. [4ji 
is (5-accurate because e and / do not overlap in x. 

A refinement is a 8 -splitting when every split point is 5-close to the curves 
that it splits (Dcf. [T2|). Although the curves in Fig. [5] are inaccurate, the refine- 
ment is a 5-splitting because e is split at points that are i5-close to it. 

Definition 12 (<5-splitting). A refinement of an approximate subdivision, S, is 
a (5-splitting if for all e £ S and e'[a', b'] £ C(e), dist(a', e), dist(6', e) < 6, where 
dist(a, e) denotes the Hausdorff distance from a point, a, to a curve, e. 

Lemma 15. The refinement algorithm output is a S-splitting (p. \20\) . 

The output error of the refinement algorithm is due to the input error and the 
refinement error. If the input is 5-accurate and the refinement is a ^-splitting, 
the output is (5-accurate. A stronger error bound is that the output has a 
realization in which each curve is <5-close to its preimage curve, called a 8- 
refinement (Def. [T3|) . A (^-refinement has a single <5-deformation that places all 
the curves in -< y order, whereas 5-accuracy allows a separate (5-dcformation for 
each pair of curves. The realizations in Fig. [7] show that the refinements in 
Fig. [5] are ^-refinements. 

Definition 13 (8- refinement). A 8- refinement of an approximate subdivision, 
S , is a refinement that has a realization such that the Hausdorff distance from 
r(e') to e is bounded by 8 for all e £ S and e' £ C(e). 

We conclude the error analysis by proving (Sec. [5]) that the stronger error 
bound is implied by the weaker one and consistency, hence that the refinement 
algorithm computes a (^-refinement of its input. 

Theorem 5. A consistent 8-splitting of a 8-accurate approximate subdivision 
is a 8-refinement (p. \20\) . 

Theorem 6. The refinement algorithm computes a 8-refinement (p. \ 21\) . 




Figure 7: Realizations of Fig. [5] that are <5-closc to the approximate subdivisions 
in Fig. |U 

3 Refinement Algorithm 

The refinement algorithm consists of four steps that enforce the four constraints 
in Def. [7] by splitting curves at endpoints of other curves. Splitting e[t, h] at p 
replaces e by a refinement, {ei[t,p], e2\Jp, h]}, such that e; -< y f (or e, y y /) iff 
e -< y f (or e y y f) and and / overlap in x. For simplicity, one can choose 
ei[t,p] and e2[p, h] to be the line segments tp and ph. Even though this may be 
inaccurate, the ultimate realization of the refinement will be close to the original 
curves (Thm.[5|). 

In the following description, a curve, e[t, h], is classified with respect to an 
endpoint coordinate, x, as follows: t x = x = h x , vertical; t x = x < h x , outgoing; 
t x < x = h x , incoming; and t x < x < h x , passing. 

Step 1 We begin with an example using curves f\ to fig and -< y equal to 
index order for curves that overlap in x (Fig. [8^). Because f-j -< y fg, w £ L(fg), 
yet because fg -< y /io, w s U(fg). Hence, curve fg and endpoint w violate 
constraint 1 if we set e = fg, p = w, and q — w. The algorithm visits, in an 
order consistent with -< y , the curves that have an endpoint whose x-coordinate 
equals w x : /i, f 2 , h, h, h, fr, fa, and f 10 . (Curve U[u,v] is vertical and 
u x = v x = w x .) It forms groups of these curves separated by passing curves: 
G\ = {/i, h, h,U, h}, passing curve / e , G 2 = {/r, fs}, passing curve fg, and 
G3 = {fio}- Each curve is added to the latest group and the constraint is 
checked. The check fails when the highest endpoint, p, in the previous group 
and the lowest endpoint, q, in the current group satisfy p y > q y . The algorithm 
splits at p each passing curve e that separate the previous group from the current 
group (because e, p, and q violate constraint 1), and it then merges these two 
groups and the split curves into a single group. In our example, the check fails 
after f w is added to G3, fg is split at w into fg and fg, G3 is merged into G2, 
and the algorithm is done (Fig. [8b) . 

The algorithm is a plane sweep with a vertical sweep line. The sweep list, 
L, consists of curves in -< y order. The events are the endpoint x-coordinates in 
increasing order. The incoming curves at x are removed from L, the non-passing 
curves are handled in -< y order as described below, and the outgoing curves are 
inserted into L. We write f -< y G or G ~< y f when / -< y g or g -< y f for every 
g in group G. The endpoints, v, of curves in G with v x = x and with minimum 
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Figure 8: Curves before (a) and after (b) constraint 1 enforcement. 

and maximum y-coordinatcs arc denoted ymin(G) and ymax(G). A curve, e, is 
handled as follows. 

1. If m = or there exists a passing / with G m < y f < y e 

increment m and set G m = {e} 
else add e to G m . 

2. While m > 1 and ymin(G m ) y < ymax(G TO _i) 3/ : 

a. For each passing curve, G rn -\ < y g -< y G m : 

i. Remove g from L. 

ii. Split g into g' and g" at ymax(G m _i). 

iii. Add g' and g" to G m _i. 

b. Insert the curves of G TO into G m _i and decrement m. 

Step 2 We begin with an example using curves f\ to f± that are incoming at 
x = bi x = b-2 x = bz x = bi x and with -< y equal to index order (Fig. |9^i). The 
pair (/2j/3) violates constraint 2 with e = fa, f = /3, b = 62 and d = 63, and 
similarly (fi,fa). We can handle (/2,/s) by splitting / 2 at 63 or by splitting 
/ 3 at 62. The latter is more accurate because dist(& 2 , /s) (shown with a dashed 
line) is less than dist(&3, / 2 ). We say that b 2 is safe for (/ 2 , /s)- Likewise, 63 is 
safe for (/i,/ 3 ) if dist(6 3 ,/i) < dist(6i,/ 3 ). We remove the two violations by 
splitting /1 into {f[,b s b±} and / 3 into {/^,& 2 &3} (Fig.^p). 

The algorithm enforces constraint 2 via safe splits of curves that violate 
it. It processes the incoming curves at each endpoint x-coordinate as follows; 
outgoing curves are processed analogously. Let s be the list of incoming curves 
in -< y order. A safe endpoint for e[w,a] G s is an endpoint, 6, of a curve, 
f[w, b] e s, with dist(6, e) < dist(a, /). All distances are measured with respect 
to the input curves: for a refinement curve e' € G(e) (Def. E]), use dist(6, e) 
in place of dist(6, e'). For endpoints a and 6, let min(a, 6) denote a ii a y < b y 
and & otherwise; likewise max(a, b). Define m(e) as the minimum endpoint over 
g >y e that is safe for e. In our example, m(/i) = &3,m(/2) = &2 7 nr(/ 3 ) = 
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(a) (b) 
Figure 9: Curves before (a) and after (b) constraint 2 enforcement. 



&3,m(/4) = 64. The algorithm computes m(e) for every e £ s then splits each 
curve, e, at the maximum m(/) over / -< y e. 

1. For each curve, e\v, a], in ,s: 

a. Set m(e) = a. 

b. Append e to /. 

c. While a y < by for the predecessor, f[w,b], of e in /: 

i. Swap e and / in /. 

ii. If a is safe for /, set m(/) = min(m(/),a). 

2. For each curve, e, in s 

Split e at the maximum of m(/) over / ^< y e in s. 

Steps 12 of the algorithm compute m(e) by inserting each curve, e[v, a], 
into an initially empty list, /, ordered by endpoint y-coordinate. First, e is 
appended to / then it is swapped with each predecessor, f[w,b], with b y > a y 
until it reaches the correct position in /. In our example, is inserted into 
I = (/i,/2) with two swaps after which / = (/3,/i,/2)- Each swap reveals a 
constraint violation, hence implies a possible update to m(/3). 

Step 3 We begin with an example with curves g[t, h], ,fi[a\, b{\ and /2[fl2, 62] 
with fi ^ y g and g ^ y fi (Fig. ITOk). Curve e = g and p = a-2 G U(e) violate 
constraint 3b with c = t, and d = b\. The algorithm constructs a monotone 
curve from t to h that is above every p £ L(g) and is below every q £ U(g). 
It visits these endpoints in increasing x order and tracks the rightmost one, 
v(g), where the curve is forced to increase or decrease. A constraint violation is 
detected when v(g) switches between increasing and decreasing, and g is split 
at the previous v(g). Specifically, v(g) starts as t, is set to 02, is unchanged at 
62 and a%, is set to 61, a violation is detected, and g is split at 02 (Fig. Hob). 

The algorithm is a sweep with L as before. Each e £ L stores an endpoint, 
v(e), and a flag, vtype(e), that indicates if v(e) is on, above, or below e. Each 
input curve, f[t, h], has a tail event at t x and a head event at h x . The events are 
handled in increasing x order with ties broken as follows: heads of non-vertical 
curves, tails of vertical curves, heads of vertical curves, and tails of non- vertical 
curves. 
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Figure 10: Curves before (a) and after (b) constraint 3 enforcement, 
tail event for f[t, h]: 

1. If / is not vertical, insert it into L with v(/) = t and vtype(/) = on. 

2. For each e € L with e -< y / and i y < v(e) y , call lower(e, i). 

3. For each g S i with g > y f and t y > v(g) y , call raise(<7,i). 

head event for /[t, /i]: 

1. For each e £ L with e ^ u / and /i y < v(e) a , call lower(e, h). 

2. For each g £ L with g y y f and /i y > ~v{g) y , call raise(<7, /i). 

3. If / is not vertical, remove it from L. 

subroutine lower(e, a): 

1. If vtype(e) = below, split e at v(e). 

2. Set v(e) = a and vtype(e) ~ above. 

subroutine raise(e, a): 

1. If vtypc(e) = above, split e at v(e). 

2. Set v(e) = a and vtype(e) = below. 

This algorithm handles strict extrema. In the general case, lower sets a 
pointer from a to v(e) if vtype(e) = above and v(e) y = a y ; likewise if vtype(e) = 
below in raise. The split steps traverse the pointers and split e at the extrema. 

Step 4 Split every vertical curve at every endpoint that violates constraint 4. 

4 Correctness 

We prove that the refinement algorithm output is a consistent refinement of the 
input (Thm. Q} and that the algorithm is inconsistency sensitive (Thm. [2]). We 
then prove that an approximate subdivision is consistent iff it has a realization 
(Thm. [3]) , which implies that the refinement algorithm output has a realization. 

Theorem 1. The refinement algorithm computes a consistent refinement. 

Proof. We prove consistency. Being a refinement follows from the definition of 
a curve split. 

Step 1 enforces the invariant that constraint 1 holds for every passing curve, 
/, and the endpoints in G\, . . . , G m . Step 2 enforces constraint 2 by assigning 
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each incoming curve an endpoint whose y-coordinate is the maximum over all 
incoming curves below it in -< y order. For step 3, let input curve e and endpoint 
q £ U(e) violate constraint 3b with c, d £ L(e); the other case is similar. If e is 
split at an endpoint, w, with c x < w x < d x before d is handled, we are done. 
Otherwise, v(e) y > c y after c is handled, so v(e) y < q y and vtype(e) = above 
after q is handled. When d is handled, vtypc(e) = above and v(e) = r for some 
r £ U(e) with r y < q y and c x < r x < d x , since no split occurs before d. Thus, 
e is split at r, which eliminates the e,q inconsistency. Step 4 trivially enforces 
constraint 4. 

Step 2 maintains constraint 1 by not changing L(e) or U(e) for a passing 
curve, e. Step 3 maintains constraint 1 by not adding an endpoint, v, to U(e), 
except by splitting e at v (proof below) . It maintains constraint 2 by splitting 
a curve, e, at the endpoint v £ L(e) or v £ U(e) with the largest or smallest v y 
value. Step 4 maintains constraints 1-3 by only splitting vertical curves. 

Proof: Suppose e -< /, v £ L(e) n L(f), and / is split at v. We show that e 
is also split at v. After v is processed, v(/) = v, since it is later split at v. By 
Lcm. [T]bclow, v(e) y < v y . But v(e) y > v y because v £ L(e), so v(e) = v. Curve 
/ is split at v when the sweep processes an endpoint, u £ U(f) with u y < v y . 
If v(e) = v then, e is split at v because u £ U{e). Otherwise, e is lowered by 
w £ U(e) n L(f) with w x < u x and e is split at v then. □ 

Lemma 1. In step 3, v(e) y < v(f) y for e, f £ L with e -< y f . 

Proof. The invariant holds initially and is preserved by sweep updates because 
if v £ L(e) overlaps / in x, then v £ L(f). If we call raise(e, v) and if v y > 
v(f) y , then we also call raise(/, v). Hence v(e) y < v(f) y remains true because 
v ( e )y = v (f)v = v - Similarly for u £ U(f). □ 

To prove inconsistency sensitivity, wc need to count constraint violations. 
A constraint 1 violation is a pair e, q such that there exists p that violates the 
constraint, and similarly the other violations are defined by the pair e, / (2), 
the pair e,p (3a), the pair e, q (3b), or the pair e,p (4). 

Theorem 2. The refinement algorithm complexity is 0(nlogn+fc) forn curves 
with k = 0(n 2 ) constraint violations. 

Proof. The definitions of the constraint violations imply that k = 0(n 2 ), since 
there are n curves and O(n) endpoints. Each curve, e, stores its poset height, 
hcight(e). For e and / that overlap in x, e -< y f iff height(e) < height (/). Each 
curve stores pointers to the curves below and above its endpoints. Splitting a 
curve updates these pointers in constant time. Sort the endpoints in x order 
with ties broken by height. This processing takes O(nlogn) time and O(n) 
space. 

Step 1 takes 0(n + k) time using a doubly linked list to represent L. Group 
operations take constant time using ymin, ymax, and the maximum curve 
height, hmax, to represent a group. 

Step 2 takes 0(n + k) time because the number of constraint 2 violations 
is O(k) after step 1, which is shown as follows. If step 1 splits e at p, p is the 
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highest point in G m _i, so the new curves satisfy constraint 2 with respect to a 
curve, / -< y e. A constraint 2 violation is created for each curve, / >~ y e, that 
has an endpoint, q, with q x = p x and q y < p y . Since q £ U(e), e and / are also 
a constraint 1 violation. 

Step 3 takes 0(nlogn + k) time. Maintain a partition of L into sublists with 
a common v and vtype. Each sublist stores pointers to its lowest and highest 
curves. Represent the partition as a balanced binary tree of sublists ordered 
by v y . By Lem. [TJ these will also be ordered by h. Store the curves in L in 
a binary tree ordered by h. To handle an event, look up the new curve in the 
latter binary tree to determine the h of its neighbors. Look up this h in the 
former binary tree to determine the sublist that will be split. The new curve 
can lower the curves below it in its sublist or can raise the curves above. Each 
case splits the sublist into lower and upper sublists. When a lower or a raise 
propagates to sublists below or above, they are merged into the newly created 
lower or upper sublists. A tail event also creates a singleton sublist for the curve 
that enters L. A head event removes its curve from its sublist and removes an 
empty sublist from the partition. 

Excluding sublist merges and curve splits, the sweep takes 0(n log n) because 
there are O(n) tree operations. Although an event can cause O(n) merges, the 
merge time is O(nlogn) because each merge decreases the number of sublists, 
whereas at most 2n sublists are created. Since each split applies to an entire 
sublist, the partition structure is unchanged and the only cost is 0(k) to update 
the curves. 

Step 4 takes 0(n + k) time, using the sorted endpoints, because the number 
of constraint 4 violations is O(fc) after steps 1-3. Each split creates at most one 
violation and is charged to a constraint 1-3 violation. □ 

We now turn to Thm. [3] We first prove necessity (Lem. [2]) . To prove suffi- 
ciency, we start with a more general sufficient condition (Def. [T51 and Lem. [3]). 

Lemma 2. If an approximate subdivision has a realization, it is consistent. 

Proof. For constraint 1, consider p £ L(e) and q £ U(e) with a x < p x = q x < b x 
(Fig. 13^). There exists a curve / -< y e with endpoint p and a curve g )^ y e with 
endpoint q. Hence, r(f) < y r(e) and r(e) < y r(g). By continuity, p < y r(e) 
or p £ r(e), but r(e) n r{f) = by Def. [5] and a x < p x < b x , hence p < y r(e). 
Similarly, r(e) < y q and hence p y < q y . 

For constraint 2, consider a x = c x (Fig. 0b). Since r(e) < y r(f), a £ r(e), 
and c £ r(f), a y < c y . The b x = d x case is similar. 

For constraint 3a, consider e[a, b] increasing (Fig. 0J;). As above, p £ L(e) 
and a x < p x < b x imply p < y r(e), and d £ U(e) and a x < d x imply r(e) < y d 
or d = b. Since r(e) is monotone and p x < d x , p y < d y . Similarly, p y < c y for e 
horizontal or decreasing. The constraint 3b proof is similar. 

Constraint 4 holds because r(e) n r(f) = by Def. [5] (Fig. @Ji). □ 
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Figure 11: (a) M(e); (b) M(f) (shaded) and r(f) (dark); (c) M'(e), r(e) (dark), 
and r(f)~ (shaded). 

Definition 14 (monotone). A set, M(e), is monotone with respect to a curve, 
e[a,b], if M(e) = e for e horizontal or vertical and otherwise M(e) is an open 
set bounded by monotone curves, 1(e) and u(e), that connect a to b (Fig. \llu ). 

Definition 15 (ordered monotone assignment). An assignment of a monotone 
set, M(e), to each curve, e, in an approximate subdivision is ordered if no M(e) 
contains a curve endpoint, and f -< y e implies M(f) < y M(e) unless e and f 
are both increasing or both decreasing in which case M(f) < y u(e) (Fig. \llb ). 

Definition 16 (p _ , p + , S~ and S + ). For a point, p, p~ = {(p x ,y)\y < Py} 
andp+ = {(p x ,y)\y>p v }. For a set, S, S~ = {J peS P~ and S+ = \J pe sP + - 

Lemma 3. If an approximate subdivision has an ordered monotone assignment, 
it has a realization. 

Proof. For each curve, e, we will choose r(e) C M(e) that will be disjoint from 
the other curve endpoints. If e is horizontal or vertical, define r(e) = M(e) = e. 
Vertical r(e) and r(f) cannot intersect because they cannot contain each other's 
endpoints. Otherwise, f < y e implies r(f) < y r(e) because M(f) < y M(e), 
except when both are increasing or both are decreasing. 

We define m for the increasing curves in an order that is consistent with < y . 
Let M'(e) equal M(e) minus r(f)~ for all increasing f < y e (Fig.lTTb). Because 
r(f) C M(f) < y u(e) (Fig.lTTb). M'(e) is connected. Since the endpoints of r(f) 
do not lie in R(e), the only part of the boundary of r(f)~ that can lie in R(e) is 
r(/), hence the lower boundary of M'(e) is monotone. Choose a monotone path, 
r(e) G M'(e), from a to 6 (Fig.QjJ). Because r(e)nr(/)~ cM'(e)nr(/)" = 0, 
r (f) <?/ r ( e ) as required. We handle the decreasing curves similarly. □ 

We complete the sufficiency proof by showing that a consistent approximate 
subdivision has an ordered monotone assignment, i?(e). 

Definition 17. For a point, p, p ++ = \p x ,oo) x \p y ,oc), p ^ = (— oo,Pa,] x 
[p y , oo), p + ~ = \p x ,oo) x (-oo, p y ], andp = (-oo,pj x (—oo,p y ]. For a set, 
A, A ++ = U P eAP ++ e ^ c - F° r t wo points, p and q, with p x < q x , A(p,q) = 

(px,qx) x (p v ,q v )- 
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Definition 18 (i?(e)). For a curve e[a,b]: 

!A(a, b) — L(eY U U (e) ^ for increasing e, 
A(a, b) — L(e)~~ U U(e) ++ for decreasing e, 
e otherwise. 

Lemmas 2] through [7] apply to a consistent approximate subdivision. 

Lemma 4. For a curve, e[a,b], if p € L(e), q <= U(e), and either p ^ {a, b} or 
q (I {a, b} 7 e increasing or horizontal implies p^ D q h = and e decreasing 
or horizontal implies p PI q ++ = 0. 

Proof. Suppose e'[a',6'] is increasing or horizontal, /' -< y e' has endpoint p' £ 
L(e'), g' >y e' has endpoint q' g U(e'), q' {a', b'}, and p' + ~ n q'~ + ^ 0, hence 
p' x < q' x and p' y > q' y . We will show a contradiction. The other cases are similar. 

If a' x < p' x = q' x < b' x , constraint 1 is violated with e[a, b] = e'[a! , b'], p = p' , 
and q = q' . If a! x = p' x = q' x , p' y < a' y by constraint 2 with e = /',/ = e', a = p', 
and c = a' . Similarly, a' y < q' , so a' y < q' < p' y < a' y and q' = a', which is a 
contradiction. Similarly, p' x = q' x = b' x is contradictory. If a' x < p' x < q' x < b' x , 
either p' y > a' y or q' y < b' y . The former contradicts constraint 3a with e = e', 
c = a, p = p' , and d = q' , and the latter contradicts constraint 3b with e = e', 
c = p, q = q' , and d = V. If a' x = p' x < q' x < b' x , p' y < a' y by constraint 2 
and constraint 3b is contradicted as before. Similarly, a' x < p' x < q' x = b' x is 
contradictory. If a' x = p' x and q' x = b' xl p' y < a' y and b' y < q' y by constraint 2, so 
b'y < q' y < p'y < a' y < b' y and q' = b', which is contradictory. □ 

Lemma 5. R(e) is monotone with respect to e. 

Proof. For increasing e, the boundary of p^ is monotone, so the boundary 
of L(eY is monotone because it is the upper envelope of a set of monotone 
curves. Similarly, the boundary of U(e) ^ is monotone. By Lem. [4j 

L(e)+- n U (e)"+ = , J peL{e)tqeU{e) p + - n q~+ = {a, b}. 

Hence, R(e) is an open set bounded below and above by the portions of the 
boundaries of L{e,y and U(e) ^ that connect a to b. Decreasing e is handled 
similarly and the other cases are trivial. □ 

Lemma 6. For e, f , and g that overlap in x, e -< y f and f -< y g implies e < y g. 

Proof. Since e and g overlap in x, e -< y g or g -< y e, but the latter implies that 
(e, e) is in the transitive closure of < y , which contradicts Def. HI □ 

Lemma 7. The assignment of R{e) to e is ordered. 

Proof. Verticals cannot contain endpoints by constraint 4. A horizontal, e[a, 6], 
cannot contain an endpoint, p, by constraint 3 with c = a and d = b. Otherwise, 
R(e) cannot contain p because one of p ,p ,P ++ has been removed. 

If e[a,6] is not decreasing and /[c, d] is not increasing and / -< y e, either 
a — c, a G U (/) — {c, d}, or c € L(e) — {a, b}. If a = c, e and / cannot both be 
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Figure 12: R(e) is shaded and R(f) has a thicker boundary. 



horizontal without one containing the other's right endpoint. We have p y < a y 
for p £ R(f) and q y > a y for q £ i?(e) with equality only for horizontal curves. 
Hence, R(f) < y R(e). If a £ U(f) - {c, d}, R(f) n a++ = 0. Since i?(e) C a++, 
<a R(e). The c £ £(e) — {a, 6} case is similar, as is e not increasing and 
/ not decreasing. 

It remains to show that if e[a, b] and f[c, d] are both increasing and / < y e, 
R(f) < y u(e), and similarly for both decreasing. It suffices to show R(f) < y q 
for each q £ U(e) for which q ^ and / overlap in x. If q x < d x (Fig. [121 q = 1), 
a curve, g ^ y e, with endpoint q also overlaps / in x, so g > y f by Lem. [6j 
q £ U(f), q ^ was subtracted from R(f), and thus -R(/) < a (7 '". If q x > d x 
(Fig. [T2j g = 2), d £ L(e) because it is an endpoint of / < y e. By Lem. SJ 

q h n d H = 0, hence d y < q y and R(f) C A(c 7 d) < y q |_. If q x = d x < b x 

(Fig. [121 g = 3), the same reasoning applies. If q x = d x = b x (Fig. [12l q = 4), 
?s/ > ^ and 6 y > by constraint 2, so R(f) C A(c, d) < y q h . □ 

Theorem 3. An approximate subdivision has a realization iff it is consistent. 
Proof. Lem. [5] proves necessity. Lems. El [SJ and [7] prove sufficiency. □ 

5 Error Analysis 

We prove that a consistent ^-splitting of a 5-accurate approximate subdivision 
is a ^-refinement (Thm. [5]) . This implies that the refinement algorithm output is 
a 8- refinement (Thm. [5]). We first prove the special case where the approximate 
subdivision is consistent, hence is its own refinement (Thm. 2]). We define the 6- 
offscts of a curve fDef. flUland Fig.[T3"k) and prove them monotone (Lem. [!)]). We 
prove that the -< y relation on the curves of a (^-accurate approximate subdivision 
implies a < y relation on their <5-offsets (Lems. [TOHTTj) . We use these tools to 
sharpen the Thm. |3] realization into one that is i5-closc to the approximate 
subdivision. 

Lemma 8. For a point, a, and an increasing (decreasing) curve, e, dist(a, e) is 
decreasing (increasing) in a x and increasing in a y for a > y e and is increasing 
(decreasing) in a x and decreasing in a y for a < y e. 

Proof. Follows from Lem. 8 in our prior paper. [7] □ 
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(a) (b) 
Figure 13: Offsets (a) and their properties (b). 



Definition 19 (e-g and e+j). The <5-offsets of a curve, e, are 

e_5 = {a < y e|dist(a, e) = 6} 
e+s = {a > y e|dist(a, e) = 5}. 

Lemma 9. The offsets e±s are monotone. 

Proof. Let a vertical line intersect e±g at a and b with a y < b y . Since dist(ct, e) = 
S and dist(6, e) = 6, dist(c, e) = S for c € [a, 6] by Lem. |SJ Hence, the line 
intersects e±s in a connected set. Likewise for a horizontal line. □ 

Lems. [10] and |TT| apply to a ^-accurate approximate subdivision. 

Lemma 10. For two curves, e and f , e < y f implies e_s < y f+g. 

Proof. By Def. [T9j e-s < y e < v e+s- For e -< y /, let a <E e_a and 6 € /+5 with 
o-x = b x (Fig. H3b'). There exist 5-deformations, e and /, with e < a /. Let a€e 
and b € f satisfy a x = b x = a x , so a y < b y . We have a y < a y by Lem. [5] because 
dist(et, e) = 5 and dist(S, e) < 5; likewise b y < b y . Thus, a y < a y < b y < b y . □ 

Lemma 11. If p <G L{e), p < y ~e+~g; if q e U{e), q > y eZg. 

Proof. We prove the p € L(e) case for e[a, b]. By construction, e < 9 e+5 hence 
{a, b} < y e+g by continuity and because dist(a, e) = dist(6, e) — < S. If e is a 
(5-dcformation of e, e = e U {a, 6} <j, e+5 and similarly e > a e_a. 

If p g L(e), p is an endpoint of / < y e and hence an endpoint of / < y e 
by Def. 111! By continuity, p y < r y for r € e with p x = r x . Since e < y e+s, 
r < y e+s hence p < y e+s- □ 

For the proof of Thm. 01 Rg(e) serves the role that R(e) did for Thm. |3] 

Definition 20 (R$(e)). For an approximate subdivision, define Rg(e) = R{e) 
for horizontal and vertical e. For increasing e, Rg(e) is R(e) minus fZg for all 
increasing f ^ y e and minus g^ s for all increasing g >i y e (Fig. \1J$ - Similarly 
for decreasing e. 
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Lemma 12. If an approximate subdivision is consistent and S- accurate, the sets 
R$(e) are monotone. 

Proof. We prove the e[a, b] increasing case. Let f[c,d] < y e be increasing. For 
q € U(f), /„5 < y q by Lem. [TTJ so f_$ < y q h , since f_g is increasing. Hence, 
f-s < y U(f) '"j which implies f_g C R(f)~, since both span the x-interval 
(c x ,d x ). Since R(f) < y u(e), R(f)~ < y u(e) and f_$ <y u(e). Likewise, 
g+s >y l( e ) f° r increasing g > y e. If / < y e, e ^j, g, and / and g overlap in x, 
f diy g by Lcm. [6] and /_£ < y g + $ by Lem. [TO] We conclude that the lower and 
upper sets that are removed from R(e) are disjoint. Furthermore, c,d i?(e) by 
Lcm. [7] and f-s <y f ■ Thus, the only part of the boundary of fZg that enters 
R{e) is part of fs, which is monotonic. Therefore, the lower boundary of Rg{e) 
is monotonic and similarly the upper boundary. □ 

Lemma 13. Let ug(e) be the upper boundary of Rs(e). If /[c, d] -< y e[a, b] are 
both increasing or both decreasing, Rg(f) < y ug(e). 

Proof. Pick p e Rs(f) with a x < p x < b x . Since Rg(f) C R(f) < y u(e), 
p < y u(e). Pick g >z y e that overlaps p in x. Thus, g overlaps / in x and g >~ y f 
by Lem. |6] Therefore, gf g was removed from Rs{f) hence p < y gj_ s . Since p 
is below the upper sets removed from Rs(e) that overlap it in x, p < y us(e). 
Hence, Rs(f) < y u s (e). □ 

Theorem 4. A consistent, 5 -accurate approximate subdivision is a 5 -refinement. 

Proof. The Rs (e) are monotone by Lem. 1121 They are an ordered assignment 
because Rs(e) C R(e), so Rs{e) < y Rs(f) for e < y f unless both are decreasing 
or both are increasing, and that case is covered by Lem. [13J By Lem. [3j there 
is a realization, r(e) C Rs{e). For p £ Rs(e), dist(p, e) < 5 by construction, so 
r(e) is a (5-dcformation of e. □ 

We prove Thm. [5] similarly, using Lem. [14] instead of Lem. [TT] and ^(e') 
(Def.HI]) instead of R s (e). 
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Lemma 14. In a 5-splitting, S' , of a S-accurate approximate subdivision, S, 
for every endpoint p £ S, curve e € S , and curve e' £ C(e), p £ L(e') implies 
p < y e+s and p £ U (e') implies p > y es- 

Proof. We prove the p £ L(e') case. Consider e'[a',b']. If p £ L(e), p < y e~^s 
by Lem. [TT] If p £ f (e), p is an endpoint of / € 5 with e -<j, /. Thus, p is 
an endpoint of a curve, /' € C(/), but by Def. El e' :< y /', so p £ U{e'). By 
Lem. HI p e {a', b'} hence dist(p, e) < S by Def. IT2l and so p < y e+5. □ 

Definition 21 (e^_ 5 , e' +(5 , i?^(e')). For increasing e'[a',b'] £ C(e), define e'_ s 
and e' +s to be e-$ and e + $ restricted to (a x ,b x ) if e is increasing; otherwise 
e'_ s = e', s = 0. Similarly for decreasing e' . Define R' s {e!) = R{e') for horizontal 
or vertical e' . For increasing e! , R' s {e') is R(e') minus f'_g for increasing 
f ^< y e' and minus g' , s for increasing g' > y e' , and similarly for decreasing e' . 

Theorem 5. A consistent 6-splitting of a S-accurate approximate subdivision 
is a 6 -refinement. 

Proof. Modify the proof of Thm.|4]by replacing Lem.llllwith Lcm.[l4]to select a 
realization r(e') C R's(e') for each curve e'[a', b'}. If e' and e are both increasing 
or both decreasing, Rg(e') lies between e-$ and e+s by construction hence r(e') 
lies within 5 of e. If e is increasing but e' decreasing, we know dist(a', e) < S 
hence a' < y e + $ hence a' + < y e+$. Similarly, b' + > y e_a. Hence, r(e') C 
a' + ~ n b'~ + lies within 8 of e. □ 

Lemma 15. The refinement algorithm output is a 5-splitting. 

Proof. Let the algorithm split a curve, e, at an endpoint, p. We show that 
dist(p, e) < 8. 

Step 1: Necessarily, p £ L(e), because p belongs to a group below e, and 
there is a violation of constraint 1. Hence, there exists q £ C/(e) with p^ = q x 
and p y > q y . By Lem. [TT1 p < y e + s an d g > y e_^. Since p y > g^, p > y as 
well. 

Step 2: For an incoming curve, e[u,a], define M(e) analogously to m(e) as 
the maximum safe endpoint over f[w,b] ^ y e. Since e ^ y e and dist(e, a) < 
dist(e,a), m(e) < y a < y M(e). Curve e is split at the maximum m(/) with 
/ < y e. Let f[w,b] achieve this maximum. Since e ^ y e, m(e) < y m(/). If 
b is safe for e (dist(6, e) < dist(/, a)), m(/) < y b from above and b < y M(e) 
by definition of M(e); otherwise, if a is safe for /, m(/) < y a by definition of 
m(e) and a < y M(e) from above. Hence m(/) < y M(e). Since both m(e) and 
M(e) are safe for e, dist(m(e), e), dist(M(e), e) < S by Lem. [T6l below. The error 
bound neglects the distance computation error, 7, which is typically negligible 
with respect to S. Since we evaluate dist approximately, we might pick the 
endpoint that is unsafe by a little bit. Strictly speaking, we should replace 5 
with S + 27. 

Step 3: We consider constraint 3 with e increasing; the other cases are 
similar. Since p £ L(e) and d £ [/(e), p < y e + s, and d > y e-s by Lem. [Til Also, 
p > y es because p x < d x , Py > d y , and es is increasing by Lem. [9] 
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Step 4: dist(p, e) = 0. □ 

Lemma 16. If e[v,a] -< y f[w,b] violate constraint 2 at a x = b x , dist(a, /) < 8 
or dist(&, e) < 5. 

Proof. Let e and / be curves that show (5-accuracy. Since e < y f and a y > b y , 
either e has a vertical segment at x = a x that contains b or / has such a segment 
that contains a. The result follows because e and / are 5-close to e and /. □ 

Theorem 6. The refinement algorithm computes a S-refinement. 

Proof. Follows from Thm. [TJ Thm. [5j and Lem. [T5] □ 

6 Set Operations 

We now turn to set operations. We compute the overlay approximate subdivi- 
sion with our sweep algorithm [7] then refine. The sweep algorithm computes 
the curve crossing points approximately and preserves the curve endpoints, so 
the output endpoints are the union of the input endpoints and the approximate 
crossing points. The overlay -< y determines the membership of the overlay faces 
in the input faces via a standard algorithm. [15] 

We modify the sweep algorithm to use red/blue sweep list insertion. When 
inserting a curve from the first input (red), find the two nearest red curves, 
using the red -< y , then compute the order of the new curve with respect to the 
intervening blue curves. Likewise for blue curve insertion. Wc avoid computa- 
tion of red/red and blue/blue orders. We also ensure that the overlay < y agrees 
with the red and blue -< y , so overlay face membership is computed correctly. 

7 Euclidean Transformations 

We transform an approximate subdivision by transforming each curve, comput- 
ing the transformed ~< y , and invoking the refinement algorithm. The computa- 
tional complexity is dominated by the refinement algorithm. 

Translating by t maps e to {p+£|p £ e} and drops each curve whose endpoints 
translate to the same point due to rounding. It does not change -< y or in- 
accuracy. Scaling (y) by w maps e to {(p x ,wp y )\p <E e} and drops the curves 
with equal endpoints. It does not change -< y , but scales ^-accuracy by if. 
Skewing by w maps e to { (p x , wp x + p y )\p <E e} and splits the skewed curves at 
their y turning points. It increases the distance between two points by at most 
a factor of z — \J w 2 + 1. Skewing a deformation yields a deformation. Hence, 
skewing scales 5- accuracy by at most z. 

Rotation Rotating by 90° maps e to {(— p y ,p x )\p G e }- The rotated curves 
are monotone because the input curves are. The rotated -< y is the horizontal 
order of the input curves. Let the rotations of e and / be e r and f r : e r < y f r 
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^x 1 ^x 

(a) (b) 

Figure 15: Horizontal cndpoint/^y for increasing (a) and decreasing (b) curves. 



iff e is left of /. It is 5-accurate because a deformation that places two curves 
in -< y order also places them in horizontal order. 

Compute the horizontal order via a horizontal line sweep whose events are 
insertion and removal of each curve at its lower and upper cndpoints. Removal 
is standard. A curve, /, is inserted by computing the order of its lower endpoint, 
v, with respect to O(logn) curves in the sweep list. Let e[t, h] be such a curve. 

1. If v x < t x , v is left of e. 

2. If / -< y e and e is increasing, v is right of e; if e is decreasing, v is left of 
e. 

3. If e -< y / and e is increasing, v is left of e; if e is decreasing, v is right of 
e. 

4. If h x < v x , v is right of e. 

Fig. [131 shows points vi, . . . ,V4 that satisfy the four rules. 

Rotating by —90° maps e to (e y ,—e x ). The rotated -< y is the opposite of 
the horizontal order of the input curves. 

Rotation by 8 is performed in two steps. Step 1 rotates the input by 90° 
and replaces 9 with - 90° for 45° < 8 < 135, rotates the input by -90° and 
replaces 8 with + 90° for —135° < 8 < —45°, and is omitted for other 8 values. 
Let s = sin.8 and c = cos 8. Step 2 scales by c, skews by s, rotates by 90°, scales 
by 1/c, skews by s/c, and rotates by —90°. The rotation formula is verified by 
direct calculation. The 5-accuracy scales by 



1 x c x \Jl + s 2 x 1 x - x \j 1 + x 1 = yi- 



This factor is bounded by y/3 because \8\ < 45° after step 1. 

8 Validation 

We have implemented set operations and Euclidean transformations for approx- 
imate subdivisions whose curves arc approximate algebraic-curve segments. A 
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(a) (b) (c) 

Figure 16: Curve representation (a), portion of b in box (b), and curve (c). 




(a) (b) (c) (d) 




(e) (f) (g) h 

Figure 17: Validation: shapes 1-4 (a-d) and their si regions (e-h). 

curve is represented by its endpoints, t and h with t x < h x , and by a monotone 
branch, b, of an algebraic curve. The curve is the portion of b in the box with 
corners t and h, plus two line segments that link t and h to it (Fig. [T6|) . If 
t y < h y and b is decreasing or h y < t y and b is increasing, the curve is the 
line segment th. When a curve is split, the new curves inherit its branch. Our 
numerical solver [7] computes points on curves by solving univariate polynomials 
and computes curve crossing points by solving pairs of bivariate polynomials. 

We validate the algorithms on sequences of set operations alternating with 
Euclidean transformations. We test the four shapes shown in Fig. [T71 The 
shapes are drawn in the box [—1, 1] x [—1,1]. The curves are line segments in 
shape 1, line and circle segments in shape 2, degree-3 curves in shape 3, and 
degree-6 curves in shape 4. We initialize so to a shape then set s^+i to the 
symmetric difference of s; with a Euclidean transform of Sj for six iterations 
(e h). The s transform is translation by (0.0123,0.0321) and rotation by 0.765 
radians. These values are halved after each iteration. 

The initial tests have high combinatorial complexity and no degeneracy. The 
program encounters no constraint violations, generates accurate outputs, and is 
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Table 1: Initial tests: n is the input size, e is the output size, S is the output er- 
ror, t is the running time in seconds, c is the constraint enforcement percentage, 
and s is the solver percentage. 



shape 


n 


e 


5 




t 


c 


s 


1 


16 


67000 


3 x 10- 


lti 


1.1 


23 


20 


2 


20 


70000 


6 x 10- 


16 


3.6 


10 


70 


3 


40 


85000 


8 x 10" 


15 


4.5 


8 


67 


4 


80 


277000 


4 x 10" 


12 


127 


1 


91 



fast. Table Q] summarizes the results. The input size, n, is the number of input 
curves. The output size, e, is the number of curves in the final, sq, subdivision. 
The output error, S, is the larger of two error metrics both of which are the 
maximum over the six iterations. The first metric is the distance between the 
monotone branch of a curve and one of its endpoints. The second metric is the 
distance between a curve, e, and an endpoint, v, with v £ L(e) and v > y e, 
or v € U(e) and v < y e. The running time, t, is for one core of an Intel 
Core 2 Duo. As the algebraic degree of the input increases from one to six, 
the constraint enforcement percentage, c, decreases rapidly and the numerical 
solver percentage, s, increases rapidly. 

We repeated the tests with the sq Euclidean transformation divided by 2 J 
for j = 1, . . . , 50. As the transformation shrinks, the six iterates of the input 
shape converge, so the symmetric differences approach degeneracy. Fig. [18] 
plots the results for the initial, j = 0, test and for the 50 subsequent tests. The 
output error and the solver percentage are omitted because they are essentially 
constant. 

The output size decreases because each iterate has fewer intersections with 
the previous iterate. The number of constraint violations is zero for j < 15, 
increases sharply from j = 25 to j = 40, then decreases to nearly zero. We 
expect that the curves become nearly identical rather than intersecting. The 
shape 3 running time increases from j = 24, peaks at double the j = time at 
j = 32, and decreases to nearly zero. The running time for the other shapes 
increases slightly if at all then decreases to nearly zero. The running time for 
shape 4 is divided by ten, so it can be displayed with the others. The percentage 
spent on constraint enforcement decreases, except for a peak in 40 < j < 50 
where the output is small relative to the number of violations. The percentage 
also decreases as the input degree increases: constraint enforcement, which is 
combinatorial, is independent of degree, whereas curve intersection, which is 
numerical, is polynomial in degree. 

9 Conclusions 

This paper presents approximate planar shape-manipulation algorithms. The 
validation results show that the algorithms are accurate and fast on both generic 
and degenerate inputs. Constraint enforcement takes a small fraction of the 



24 



e k 




Figure 18: All tests: j is the iteration, e is the output size, k is the number of 
constraint violations, t is the running time, and c is the constraint enforcement 
percentage. 
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running time, whereas numerical computation takes most of the time. The 
results support our thesis that inconsistency sensitive algorithms are efficient 
and accurate. 
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